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Abstract 

This paper is devoted to study the possibility of forming anisotropic 
compact stars in modified Gauss-Bonnet, namely called as f(G) the¬ 
ory of gravity which is one of the strong candidates, responsible for 
the accelerated expansion of the universe. For this purpose, we have 
used analytical solution of Krori and Barua metric to the Einstein 
field equations with anisotropic form of matter and power law model 
of f{G) gravity. To determine the unknown constants in Krori and 
Barua metric, we have used the sample of compact stars, 4171820-30, 

Her X-l, SAX J 1808-3658. The physical behavior of these stars have 
been analyzed with the observational data. In this setting, we have 
checked all the regularity conditions and stability of the compact stars 
4171820-30, Her X-l, SAX J 1808-3658. 
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1 Introduction 

General relativity (GR) describes many interesting phenomena of gravity as 
weak force using the classical gauge theory. Although it is well stablished 
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both for gravity and to give us an ad hoc model of Universe, there are two 
regimes of physics which we need to modify or motivate new model. The first 
area is cosmology, the dynamics of the universe in large scale, where obser¬ 
vational data indicates that we live in an accelerating universe (Riess 1998, 
2000 ) and (Perlmutter 1998). There is no simple and fully satisfying way 
to address this acceleration using the formalism of GR. If we keep GR, we 
need to address it to an exotic matter. The only possible and reasonable way 
is to modify GR by geometrical corrections, is so called as modified gravity 
(Nojiri 2010). Another regime is the ultraviolet, when we need to modify 
graviton’s propagator appropriately, specially when we need to present a 
model for quantum gravity. In this case, we need to handle a new formalism 
to cover Lorentz non-invariance behavior. A way is modification of GR in 
an adequate form by adding higher order corrections to the Einstein-Hilbert 
action. Different types of modifications for GR have been proposed in litera¬ 
ture: the Erst and the simplest one was obtained by replacing the Ricci scalar 
R by an arbitrary function f(R), this modification is one of the oldest one 
and originally proposed in (Buchdahl 1970) and later revisited by Jedamzik 
et al. (2010) and Martin et al. (2014) to address the cosmic acceleration 
and also to solve some problems on early universe, specially inflation. It has 
been demonstrated (Capozzicllo and De Laurentis 2011, Nojiri and Odintsov 
2011) that f(R ) gravity works as a realistic and physically trustable model. 

Curvature corrections can be included in the other forms. For example 
one can consider a non minimally model of gravity in the form of f(R, L m ) 
where L m denotes the matter Lagrangian and R is the scalar curvature 
(Harko 2010). Another modification is f(R , T/) where T^ u is energy-momentum 
tensor of matter fields. There are several applications of the above mentioned 
models in literature. From string theory point of view, the next higher order 
correction to the GR can be induced by Gauss-Bonnet topological invariant 
is defined by the following: 

r< _ c>2 a t) tdhv , 75 ryfiuXa 

Ct — It — 4 + rC^Xo-lt 

Inspired from high energy physics, by looking on GR as the low energy limit 
of an affective quantum theory, string theory a valid modification to the 
GR was poropsed in which we replace R by R + /(G) (Nojiri et al. 2005). 
As we know GB term in four dimensions has no dynamics so a way is to 
couple it with the matter held or to add it by a non linear form /(G). 
Cosmology of early and late time universe has been studied as well by GB 
term (Capozzicllo 2014 and Nojiri 2006). Specially exact solutions to this 
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GB-moclified gravity (Nojiri 2000, 2001, 2002) evaporation (Sebastiani et 1. 
2013) and cosmic strings (Houndjo et al. 2013 and Rodrigues et al. 2014) 
have been investigated in literature. 

A class of exact solutions of Einstein field equations with anisotropic 
source was presented by Mak and Harko (2004). These solutions show that 
energy density as well as tangential and radial pressure are finite and positive 
inside the anisotropic star. Chaisi and Maharaj (2005) formulated an algo¬ 
rithm for the solutions of the field equations with anisotropic matter source. 
By using Chaplygin gas equation of state (EOS), Rahaman et al. (2012) 
extended the Krori-Barua (1975) solution to the charge anisotropic source. 
Recently (Kalarn et al. 2012, Kalam et al. 2013) there has been growing in¬ 
terest to study the models of the compact objects using the anisotropic source 
with Krori-Barua metric. Hossein et al. (2012) formulated the anisotropic 
star model with varying cosmological constant. Bhar et al. (2015) discuss the 
possibility for the existence of higher dimensional compact star. Maurya et 
al.(2015) formulated the charged anisotropic solutions for the compact stars. 

Compact stars are relativistic massive objects which can be described by 
GR as well as modified gravity (Camenzind 2007). Due to the small size and 
extremely massive structure they have strong gravitational force. Recently, 
there has been growing interest to study the compact stars in modified theo¬ 
ries of gravity like f(R) and f(G ) (Astashenok et al. 2014a, 2014b). Specif¬ 
ically, in the last reference (Astashenok et al. 2014a), the authors presented 
neutron stars solutions for viable models of f(G) gravity. In this paper , mo¬ 
tivated by the works of (Morneni et al. 2014) and (Astashenok et al. 2014a), 
we study compact stars solutions and their dynamical stability for a viable 
power law model of f(G). We will construct exact solutions for stars which 
are comparable with observational data. Our plan in this work is as follows: 
In Sec.2, we present basic equations of motion for f(G) gravity. The analytic 
solution for the viable f(G) model is presented in Sec.3. Sec.4 deals with the 
physical analysis of the given system. The last section summaries the results 
of the paper. 
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2 Equations of motion of compact stars in 
modified Gauss-Bonnet gravity f(G) 


To study the dynamics of compact stars in modified GB gravity /(G), we 
need to specify our gravitational model. The suitable form of GB action was 
presented as the following (Nojiri et ah 2010): 
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( 1 ) 


Where , G is the GB term, R Ricci scalar, k 2 = 8n( ^ N denotes gravitational 
coupling constant, Gjv the usual Newtonian constant, c velocity of light in 
vacuum space time, S m is matter action. We assume that our GB action 
given by ([1]) is defined for a viable form of /(G) which is consistent with 
the observational data in accelerating universe,from different observational 
constraints, like solar system tests, Cassini experiment and so on. Further¬ 
more, we assume that /(G) is smooth function of argument G and have all 
higher order derivatives / n - 2 (G). Using the matter action S m , it is possi¬ 
ble to define energy-momentum of matter fields by the standard definition 
T /y = — • If we take metric g pu as dynamical variable, the full set 

of the equations of motion derived from (JT]) are written in following form 
(Nojiri et al. 2010): 


R/j,u ppvcr RpisQcrp R'paQvp, RpvQcrp “ 1 “ RpaQvp 

+ Y (STa'Wp - 9 iut 9 up ) V'V'/g + (G/g - /) g, lv = K 2 T jlv 


( 2 ) 

( 3 ) 


here fee... — d /g^ ■ We use the convention for curvature tensor by adopting 

the signature of the metric g pu as (H-) , and consequently covariant 

derivative for a vector field is defined by V/W = d p V u — T X U V\ and the 
Riemann tensor is given by R a pvp = d u T° p - d p T + T pp T° u - F^T^ . 
Matter sector satisfies an extra conservation law V ll T /w = 0. 

To have proceed on stars, we choose the metric of the compact star as 
static (time independent) and non rotating , spherically symmetric in normal 
coordinates = ( ct , 0, <p), as the following form : 


ds 2 = c 2 e 2<t> dt 2 - e 2X dr 2 - r 2 (dd 2 + sin 2 6d(p 2 ). (4) 
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We assume that matter sector is filled by a fluid with the non zero com¬ 
ponents of the anisotropic energy momentum tensor are given by T” = 
diag(pc 2 , —p r , —pt, ~Pt ) where p is energy density, p r radial pressure and 
Pt tangential component of the pressure. If we write down the equations of 
motion given in (|3]h we obtain the following set of independent equations for 
(/i, v) = (ct, ct ) and (p, v) = (r, r): 

J_(2rA' 2A - 1) + 8e- 2A (/ GG (G" - 2A'G') + / GG g(G" 2 ) ~ 2(0" + 0' 2 )) 

+ (Gf G ~ f)e 2X = k 2 pc 2 e 2X T (5) 

J_(2n^ + 1) - (Gf G - f)e 2X = K 2 p r e 2X . (6) 

Trace of (J3J) gives us the following auxiliary equation: 

R + 8 G pa V p V"f G - 4 (Gf G - /) = -n 2 (pc 2 - p r - 2/y). (7) 

If we expand it using the metric flU, we obtain: 

<t>" + ^ - </>'X + -«>'- A') + W/fi + 8e- 2A + W — 

x (Jgg(G" — 2A , G") + /ggg(G / 2 j + 4(G/ g — f)e 2X = n 2 e 2X (pc 2 - 


Pr - 2/?*). 
( 8 ) 


If we specify the form of the metric functions, then we are able to calculate 
pressures and density using Eqs.(l5i)-(lg|). 


3 Solutions for a viable f(G ) model 

In this section, we shall try to solve (I5ll6ll8l) for a given viable model of f{G). 
For this reason, we consider the power law model /(G) = aG n proposed in 
(Cognola et ah 2006), where a is arbitrary constant and n > 0. If n < |, the 
/(G) term dominates on the Einstein term in the regime of small curvature. 
We parameterize the metric as the following (Krori and Barua. 1975): 

2A = Ar 2 , 2/ = Br 2 + G, 
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where A, B and C are arbitrary constants to be evaluated by using some 
physical matching conditions (see Table 1). Indeed this solution was treated 
as a singularity free exact solution for astatic charged fluid sphere. So our 
main motivation to introduce this set of functions for static metric goes back 
to the singularity free structure of our extended model for compact star. It 
is easy to show that neither densities nor curvature terms become singular 
in the vicinity of r — 0 with this appropriate set of metric functions. 

Energy density and pressure profiles are given by the following expres¬ 
sions: 


,-Ar 2 T 


Pr 


K 


-A;(ir 2B - e Ar2 + 1) + aG n (n - 1) 


( 9 ) 


n —Ar 2 


P 


Asz> - s 

(2r 2 A + e Ar2 - 1) + 2-— ( an(n - 1 )G n ~ 3 (G" - 2ArG') (10) 

K 2 \ 


-2 Ar 2 


n 2 r 2 


+an(n — l)(n — 2)G n 3 {G') 2 


1 — e 


Ar 2 


a 


2B(l-r 2 ) ) +2— G n {n-{M 


Pt 


e ^ '~„2 A , „Ar 2 - 4e 


-2Ar 2 


2 K 2 T 2 


(2 r 2 A + e Ar " - 1) + 


K 


(an{n - l)G n - 3 (G'" + Ard) (12) 


+an(n — l)(n — 2)G n 3 (G) 


1 — e 


Ar 2 


2 (B - BA) | - 2AG"(n - 1) 


,-Ar 2 


(B + B 2 r 2 — ABr 2 H —(Br — Ar) + 


l-e / 


4e 


-2 Ar 2 


9 V 1 J 1 9 

k r r 




2 Br 


1 — e 


Ar 2 


Ar 2 


2k 2 


) (an(n. - 1 )G n ~ 1 (G" - 2Ard) + an(n - 1 )(n - 1 )G n ~ 2 {G') 2 ) 

(13) 


-4(4 r 2B - e Ar2 + 1) + aG n (n - 1) 


Also, the radial and transverse equations of state are dehned by the following 


6 



equations: 


W r 


Pr . e 
p 


—Ar 2 


X 


K 
Ar 2 


--^(4r 2B - e Ar2 + 1) + aG n (n - 1) 


(14) 


{2r z A + e Ar2 - 1)+ 


n 2 r 2 


8e 


-2 Ar 2 


+an(n — 1 ){n — 2)G n A (G ) 


K 2 

1 — e 


(cm(n - l)G n_3 (G" - 2ArG / ) 


Ar 2 


a 


- 2(S - Br z ) + — G n (n - 1) 


n- 


-l 


(15) 


w t = 


e Ar2 o . Ar 2 -' 4e 


-2Ar 2 


2 K 2 r 2 


(2r 2 A + e Ar - 1) + 


K 


(an(n - 1 )G“- 3 (G" + .44,'') 


(16) 


+an(n — l)(n — 2)G n A (G ) 


—Ar 2 o 

e ' - -i2J 4 d „2 z 


1 — e 


Ar 2 


- 2(B - Br 2 ) | - |4G”(» - 1) 


+ 


k 2 

1 - e Ar2 


( B + Err 2 — ABr 2 H —(Br — Ar) + 


2 K 2 

/ 

l-e Ar \ 4e- 2Ar2 /2Sr 


/r 


n —Ar 2 


cm(n - l)G n (G* - 2 ArG ) + an(n - 1 )(n - l)G n ' 2 (G ) 


+ 


x 


2k 2 


- 2 - 


A- (4r 2B - e Ar ' + 1) + aG n (n - 1) 


Ar 2 


zip—/ 

(2r 2 7l + e Ar2 - 1) + 2-— ( an(n - 1 )G n ~ 3 (G" - 2ArG!) 

K 2 \ 


-2 Ar 2 


2 k 2 t 2 

+an(n — l)(n — 2)G n ~ 3 (G) 2 


1 — e 


Ar 2 


a 


- 2 (B - Br 2 ) + 2 — -G n (n - 1) 
2k 2 


-i 


(17) 


The behavior of density, radial and transverse pressures and equation of state 
(EoS) parameters are given in figures 1 - 5 . In this case EoS parameters have 
values 0 < u t < 1 and 0 < ay < 1, which shows the fact that star consists 
of ordinary matter and effect of f(G) model in the present setup. We would 
like to mention that for the plotting of graphs the values of A and B have 
been taken from the table 1. 
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Figure 1: First, second and third graphs represent the density variation of 
Strange star candidate Her X-l, SAX J 1808.4-3658(SS1) and 4U 1820 - 30, 
respectively. 
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Figure 2: First, second and third graphs represent the radial pressure vari¬ 
ation of Strange star candidate Her X-l, SAX J 1808.4-3658(SS1) and 4U 
1820 - 30, respectively. 
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Figure 3: First, second and third graphs represent the transverse pressure 
variation of Strange star candidate Her X-l, SAX J 1808.4-3658(SS1) and 
4U 1820 - 30, respectively. 
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Figure 4: First, second and third graphs represent the EOS parameter uj r 
variation of Strange star candidate Her X-l, SAX J 1808.4-3658(SS1) and 
4U 1820 - 30, respectively. 
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Figure 5: First, second and third graphs represent the EOS parameter uj t 
variation of Strange star candidate Her X-l, SAX J 1808.4-3658(SS1) and 
4U 1820 - 30, respectively. 
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4 Physical Analysis 

In this section, we shall discuss following features of our model: 


4.1 Anisotropic Behavior 

From Eqs. (J5]) and dHJ), we get the radial gradient of pressures: 


dp r 

dr 


1 
k 2 . 


2e~ Ar2 ( 1 - e Ar2 


+ Ar 2 ^1 + (AB + aG n (n — l))jj (18) 


= ^2 ( 2e ' 3Ar2 ( e 2Ar2 ( - 1 + e Ar2 + 2Ar^j + Ae 2Ar2 r 2 (e Ar2 2 Ar - l) - (19) 

4e 2Ar2 (2 + e Ar2 ) + 16 AanG n ~ 3 (n - l)r 2 (l + e Ar2 ( - 1 + 2 (B - Sr 2 )r 2 ) G' 

(—2 + n — A') + G") - 8 anG n ~\n - 1) (e Ar2 - Ar 2 - l) (g'( n - 2 - A') + G") 
-4cmG n “ 3 (n - l)r(l + e Ar2 (- 1 -2(/i- Br 2 )r 2 )) ((n - 2 - A')G"A" + G®))) 
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We compute second order derivatives as the following: 


% =^[e- 3 * a (-2e 2 * J (3(-l + e* 1 ) + ^ 2 (-3-2^ 2 (20) 

+2v4r 2 (-1 + 2Ar 2 ))) - 8 anG n ~ 3 (n - 1) (W' + G' ^6(1 + 3Ar 2 (l + 2 Ar 2 )') 

^ — 2 + n — X'^j + 4r(l + 3Ar 2 )\" — r 2 A (3) + e Ar2 (2(n — 2) ^ — 3 + 2 Ar 2 (^ — 3 
+2 (—B + Br 2 )r 2 + AAr 2 (^ — 1 + 2 (B — 5r 2 )r 2 ^j j + ^6 + AAr 2 (?> + 2 (B — Br 2 )r 2 
+4Ar 2 ^l + 2 (—5 + 5?’ 2 )r 2 j j j + 4r^ — 1 + 2Ar 2 (^ — 1 + 2 (B — 5r 2 )r 2 j j X" 

+r 2 ^1 + 2 (-B + Br 2 )r 2 ^j A®^ + ri^ZG" — 2n + 3Ar(3 + 4r — 2 nr 
+6 Ar 2 ^ + (2 + 6Ar 2 )X' — rX" j — ^4 + 2r — nr + 12ALr 2 + rX'^j G® + rG^ 1 ^ 

+e Ar2 (^2G" ^ — 3 + 2 (n — 2)r + 2ALr 2 (^ —3 + 2r( —2 + n + ( —B + Br 2 )r 
—2 (B — Br 2 )(n — 2)r 2 j + AAr 2 (^ — 1 + 2 (B — Br 2 )r 2 ^ + r(^ — 2 + AAr 2 (^ — 1 
+2 (B — Br 2 )r 2 ^ X' + r^l + 2(— B + Br 2 )r 2 A"^ + r(^4 + (n — 2)r^ — 1 
+2 (B - Br 2 )r 2 ^j + 3Ar 2 (l + 2 (-B + Br 2 )r 2 ) + r(l + 2(-5 + 5r 2 )r 2 ) A')g ( 3) 
+r(-l + 2^S- J Br 2 )r 2 )G (4} )))) (21) 


^ [2e- Ar2 (3 ( - 1 + e Ar2 ) + Ar 2 ( - 3 + (45 + aG n (n - 1)^) 

-2Ar 2 (l + (45 + «G> - l))^ 2 ))) (23) 


We observe that at center r = 0, our model provides that 


dp 

= 0, 

dp r 

dr 

dr 

d 2 p 

< 0, 

d 2 p r 

dr 2 

dr 2 


(24) 

(25) 


This indicate maximality of radial pressure and density. This fact implies 
that p and p r are decreasing function of r as shown in figures 1,2 for a class 
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of strange stars. The measure of anisotropy is 


A = -(pt-pr), 

r 


which takes the form 
A 


2 / 2 r 

= -{Pt-Pr) = -5- 


n 2 r 


-,—Ar 2 


~ &A + S--D+ k2 


4e- 2 ^ 2 


(. an(n - 1 )G n ~ 3 (G" + \'g') 


+an(n — 1 )(n — 2)G n 3 (G' , ) 2 j 


1 - e y 


a 


2 (B - Br 2 ) ) - 3—G> - 1) 


_^^2 q -j ^^2 

e —— (B + S 2 r 2 - ABr 2 + -(Br - Ar) + —|—' 

ry ry A 


4e 


-2Ar 2 


A: 2 


2£>r 


1 — e 


Ur 2 


|an(n — l)G n 1 (G ,,/ — 2ArG') + an(n — l)(n — 1 )G n 2 (G ,, ) 2 j 


D -Ar 2 


— Ar 2 


2k 2 


( - yS Ar2B - eAr 2 + !) + - !)) - ( 2 ^y- ( - ^( 4r2B - ^ + !) 


+«G"(n - 1) 


(26) 


The anisotropy will be directed outward when p t > p r this implies that 
A > 0 and directed inward when p t < p r implying A > 0. In this case 
A > 0, for larger values of r for a class of strange stars as shown in figures 
6. This implies that anisotropic force allows the construction of more mas¬ 
sive stars. In order to comprehend some general results associated with the 
strong gravitational fields, we include weak energy condition (WEC), null en¬ 
ergy condition (NEC), strong energy condition (SEC) and dominant energy 
condition (DEC). For an anisotropic fluid, these are defined as 


NEC: p + p r > 0, p + pt> 0, 

WEC: p > 0, p + p r > 0, p + Pt > 0, 

SEC: p + p r > 0, P + Pt> 0, p + p r + 2pt>0, 

DEC: p>\p r I, P > \pt\- 


We find that our model satisfies these conditions (see Figs.6-8)for specific 
values of mass and radius which helps to find the unknown parameters for 
different strange stars. 
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Figure 6: Evolution of energy constraints at the stellar interior of strange 
star Her X-l. 
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Figure 7: Evolution of energy constraints at the stellar interior of strange 
star SAX J 1808.4-3658. 
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Figure 8: Evolution of energy constraints at the stellar interior of strange 
star 4U 1820-30. 


4.2 Matching Conditions 

The intrinsic metric of the boundary surface will be the same whether it is 
derived from the internal or external geometry of the star. It confirms that 
for any coordinate system the metric tensor components will be continuous 
across the boundary surface, for this purpose the matching conditions are 
required for the interior metric (j4j) , so that the system of field equations can 
be solved constrained to the boundary conditions that radial pressure p r = 0 
at r = R. The interior metric ([4]) can be matched smoothly at r = i?, to an 
exterior Schwarzschild metric given by 

d S 2 = ^ ^ dt 2 _ ^ ^ 1 dr 2 _ r 2^ 0 2 + gin 2 (37) 

At the boundary r = R continuity of the metric functions gtt , g r r and at 
the boundary surface yield, 


9tt 9tt 5 


dg tt _ dgtt 

dr dr ’ 


(28) 
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where — and +, correspond to interior and exterior solutions. From the 
interior and exterior matrices, we get 


A 

B 

C 



(29) 

(30) 

(31) 


For the given values of M and R for given star, the constants A and B are 
given in the table 1. We would like to mention that the values of M and R 
have been taken from Lattimer et ah (2014), Li et ah (1999), Leahy et ah 
(2007) to calculate the values of A and B. 


Table 1: Values of constants A and B for the given values of M and R of 
compact star 


Strange Quark Star 

M 

R(kmn) 

MM 

A(km 2 ) 

B(km 2 ) 

Her X-l 

O.88Af 0 

7.7 

0.168 

0.006906276428 

0.004267364618 

SAX J 1808.4-3658 

1.435M 0 

7.07 

0.299 

0.01823156974 

0.01488011569 

4U 1820-30 

2.25Af 0 

10.0 

0.332 

0.01090644119 

0.009880952381 
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4.3 Stability 

Here, we discuss the stability of strange stars for this purpose the sound 
speeds V 2 v and V 2 v along along radial and transverse directions are given by 

^ f s ~ *** +Ar2 ( 1+ < 4B +“ G ” ( " -!»))] x hh 

x (2e~ 3Ar2 (e 2Ar2 ( - 1 + e Ar2 + 2Hr 2 ) + He 2Ar V 2 [e Ar °~2Ar - l) 

-He 2Ar2 (2 + e Ar2 ) + 16HcmG n_3 (n - l)r 2 (l + e Ar2 ( - 1 + 2(B - Br 2 )r 2 ^j G' 

(-2 + n - A) + G") - 8 anG n ~ 3 (n - 1) (e Ar2 - Ar 2 - l) (g'( n - 2 - A) + G") 

-4cmG n ' 3 (n - l)r(l + e Ar2 ( - 1 + 2(5 - Br 2 )r 2 )) ((n - 2 - A)G"A + G (3) ))) _1 , 

(32) 

x ^2e _3Ar2 (e 2Ar2 ( - 1 + e Ar2 + 2Hr 2 ) + He 2Ar V 2 (e Ar2 2Ar - l) 


-4e 2Ar2 (2 + e Ar2 ) + 16AanG n ~ 3 (n - l)r 2 (l + e Ar2 ( - 1 + 2 (B - Br 2 )r 2 ) G' 

(—2 + n — + G") — 8 anG n ~ 3 (n — 1) ^e Ar2 — H?’ 2 — lj [Cj (n — 2 — A) + G"^j 

-4 anG n ~ 3 (n - l)r(l + e Ar2 ( - 1 + 2 (B - Br 2 )r 2 )) ((n - 2 - A 2 )G"A + G (3) ))) 


From above equations, we have 
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V 2 - V 2 v = [[ [c 2 e 2Ar2 ( - 2 (l - e Ar2 (l + [AB + aG”) (n - 1)) r 2 ) ) + e Ar2 

x ( - 1 + e" Ar2 - 2 Ae~ Ar2 r 2 + e~ Ar2 ( - 1 + e Ar2 + 2Ar 2 ) + Ae“ Ar V 
( - 1 + e Ar2 + 2Ar 2 ) - 2Ae~ 2Ar2 r 2 ( - 1 + e Ar2 (l + r ( - 2Ar + 2 Br - ABr 3 
Br + 5 2 r 3 ) ) )) - 2e" 2Ar2 (l + e Ar2 ( - 1 + Ar 2 - Br 2 + 5 2 r 4 ) + Ar 2 ) 

-8 ane- 2Ar2 G n ~ 3 {n - 1) (l + e Ar2 (-l + Ar 2 )) ((n - 2) (g')* - G'A - G") 
+16 Aane~ 3Ar2 G n - 2 (n - l)r 2 ( - 1 + e Ar2 (l + 25r 2 )) (2GarG' - (n - 1 )(G') 2 
-GG") + 2 ane- 2Ar2 G n ~ 3 (n - 1 )r(aG n (n - l)r 2 - 2k 2 ( - 1 + e Ar2 
+2 - GG") ) + (5 - Br 2 )r 2 ( - AG" + G' ( 2(71 - 2)G" - a) G (3) ) 
x 4cme _3Ar2 G n_2 (n - l)r(- 1 + e Ar2 (l + 25r 2 )) (g(G ( 3) - 2ArG') 

+2G' ((n - 1)G" - GArG') ) 

x (2(e 2Ar2 (—1 + e Ar2 + 2Ar 2 ) + Ae 2Ar V 2 (-l + e Ar2 + 2Ar 2 ) - Ae 2Ar2 r 2 
(2 + e Ar2 ) + 16AanG”- 3 (n - l)r 2 (l + e Ar2 (-l + (5 - 5r 2 )r 2 )) 

(G'(—2 + n - A) + G") - 8 anG n ~ 3 (n - 1) ( - 1 + e Ar2 - Ar 2 ) 

(G\-2 + n - A)G") - 4 anG n ~ 3 (n - l)r(l + e Ar2 ( - 1 + 2 
(5-5r 2 )r 2 ))^ —2 + n —A)G" -G'A + G (s) ))) 


From figure 6, we can see that | v 2 t — v 2 r |< 1. This is used to check whether 
local anisotropic matter distribution is stable or not. For purpose, we use 
the cracking concept introduced by Herrera (1992) which can examine that 
potentially stable region is such region where radial speed of sound is greater 
than the transverse speed of sound. Hence, our proposed compact star model 
is stable. 
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Figure 9: First and second graphs represent the variation of A and V 2 v — V r 2 v 
respectively. 
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5 Conclusion 


On the basis of cosmological observations it has been concluded that that 
there are two phases of accelerated expansion in our universe: cosmic inflation 
in early universe and acceleration in the current expansion of the Universe. 
The investigation of current cosmic expansion and nature of DE has been 
widespread among the scientists. For this purpose, numerous efforts have 
been made based upon different strategies to modify the General Relativity. 
The /(G) gravity is one of the modifications of General Relativity. 

The present paper deals with the formulation of analytical models of 
compact stars with the anisotropic gravitating static source in the frame¬ 
work of /(G) gravity. For this purpose, we have considered the power law 
form of /(G) gravity model, further the stars are assumed as anisotropic in 
their internal structure. The analytic solution in /(G) gravity have found by 
matching the interior metric with the well-known exterior metric. The phys¬ 
ical analysis of the results implies the following properties of the anisotropic 
compact stars in /(G) gravity: 

• The bound on the EoS parameters are given by 0 < uy < 1 and 0 < 
l o r < 1, which is consistent with normal matter distribution in /(G) 
gravity. 

• The density and pressures attain the maximum value at the center. 
Hence, these are decreasing functions. 

• It has been found that the anisotropy will be directed outward when 
Pt > p r this implies that A > 0 and directed inward when P t < P r 
implying A < 0. In this case A > 0, for larger values of r for various 
strange stars. This implies that anisotropic force favors the construc¬ 
tion of more massive star in /(G) gravity. 

• The variation of V 2 v — V 2 v of different strange stars is shown in Figure 
6 . This implies that difference of the two sound speeds, i.e., V 2 v — V r 2 v 
satisfies the inequality \V^ — V 2 v \ < 1. Thus, our proposed strange star 
model is stable in /(G) gravity. 
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Response of the Referee’s Report 
about the paper entitled 
Anisotropic Compact Stars in /(G) 

Gravity 

We have inserted the relevant explanation in the revised version in the 
light of referee’s comments. These are given in the following: 

(Response-Major Corrections) 


1. All the indicated and many other typos have been removed 

2. Here, we have provided the references i.e., and later revisited by Jedamzik 
et al. (2010) and Martin et al. (2014). 

3. Here we have also provided reference i.e., It has been demonstrated 
(Capozzicllo and De Lanrentis 2011, Nojiri and Odintsov 2011) 

4. For the discussion of Anisotropic compact stars we have presented the 
following paragraph at page 3 and have mentioned the required refer¬ 
ence 

”A class of exact solutions of Einstein field equations with anisotropic 
source was presented by Mak and Harko (2004). These solutions show 
that energy density as well as tangential and radial pressure are finite 
and positive inside the anisotropic star. Chaisi and Maharaj (2005) 
formulated an algorithm for the solutions of the field equations with 
anisotropic matter source. By using Chaplygin gas equation of state 
(EOS), Rahaman et al. (2012) extended the Krori-Barua (1975) so¬ 
lution to the charge anisotropic source. Recently (Kalarn et al. 2012, 



Kalam et al. 2013) there has been growing interest to study the models 
of the compact objects using the anisotropic source with Krori-Barua 
metric. Hossein et al. (2012) formulated the anisotropic star model 
with varying cosmological constant. Bhar et al. (2015) discuss the 
possibility for the existence of higher dimensional compact star. Mau- 
rya et al.(2015) formulated the charged anisotropic solutions for the 
compact stars.” 

5. For the Krori-Barua ansatz at page 5, we have added the following 
motivation 

” Indeed this solution was treated as a singularity free exact solution for 
astatic charged fluid sphere. So our main motivation to introduce this 
set of functions for static metric goes back to the singularity free struc¬ 
ture of our extended model for compact star. It is easy to show that 
neither densities nor curvature terms become singular in the vicinity of 
r = 0 with this appropriate set of metric functions.” 

also, for the ansatz given by Eq.(27) at page 13, we have added 

The intrinsic metric of the boundary surface will be the same whether 
it is derived from the internal or external geometry of the star. It 
confirms that for any coordinate system the metric tensor components 
will be continuous across the boundary surface, for this purpose the 
matching conditions are required for the interior metric (|4]), so that the 
system of field equations can be solved constrained to the boundary 
conditions that radial pressure p r = 0 at r = R. The interior metric 
(2]) can be matched smoothly at r — R, to an exterior Schwarzschild 
metric. 

6. All the plots in the Figs. 1, 2, 3, 4, 5 and 9 have been enlarged for 
clarity 

7. Energy conditions for the different compact stars have been shown in 
Figs. 6-8. 

8. We have modified the table 1 by adding few column, data set from 
Lattimer et al. (2014) have been used for masses and radii of the stars. 
The values of the constants in table have been calculated from our 
models. The references are only for observed values of mass and radius 
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of each star. To, avoid any confusion it all has been clarified before 
table 1. 

Finally, we would like to thank the Referee for his keen interest to study 
the paper in detail. Also, his fruitful comments are highly acknowledged. 
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